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Abstract 

In this Letter, by reconstructing the Om diagnostic and the deceleration parameter q from 
the latest Union2 Type la supernova sample with and without the systematic error along with 
the baryon acoustic oscillation (BAO) and the cosmic microwave background (CMB), we study 
the cosmic expanding history, using the Chevallier-Polarski-Linder (CPL) parametrization. We 
obtain that Union2+BAO favor an expansion with a decreasing of the acceleration at z < 0.3. 
However, once the CMB data is added in the analysis, the cosmic acceleration is found to be 
still increasing, indicating a tension between low redshift data and high redshift one. In order 
to reduce this tension significantly, two different methods are considered and thus two different 
subsamples of Union2 are selected. We then find that two different subsamples+BAO+CMB 
give completely different results on the cosmic expanding history when the systematic error 
is ignored, with one suggesting a decreasing cosmic acceleration, the other just the opposite, 
although both of them alone with BAO support that the cosmic acceleration is slowing down. 
However, once the systematic error is considered, two different subsamples of Union2 along with 
BAO and CMB all favor an increasing of the present cosmic acceleration. Therefore a clear-cut 
answer on whether the cosmic acceleration is slowing down calls for more consistent data and 
more reliable methods to analyze them. 

PACS numbers: 95.36. +x, 04.50.Kd, 98.80.-k 
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I. INTRODUCTION 



The fact that our Universe has entered a state of accelerating expansion at redshifts 

Dservational data, including 



3], the cosmic microwave 



less than ~ 0.5 is well established by various independent o 
the Type la Supernova (SNIa) [lj], the large scale structure 
background (CMB) radiation |4j, and so on. In order to explain this observed phenomena, 
one usually assumes that there exists, in our Universe, an exotic energy component, named 
dark energy (DE), which has negative pressure and thus can generate a repulsive force. 
It dominates our Universe and drives it to an accelerating expansion at recent times. 
Since the equation of state (EOS) w of dark energy embodies its properties, one may 
adopt a parametrized form of w(z) with several free parameters, such as the Chevallier- 
Polarski-Linder (CPL) parametrization [5], to probe the cosmic expanding history and 
the evolutionary behavior of dark energy from observations. 

However, the results are different, sometimes even contradictory, when different ob- 

12]. For example, by investigating the diagnostic Om jg], 



servational data are used 
which is defined as 



10 



Om(z) 



I 



1 + z) 



E{z) = H(z)/H , 



(1) 



and the deceleration parameter q from the Constitution SNIa [7] along with the baryonic 
acoustic oscillation (BAO) distance ratio data js), and using the CPL parametriza- 
tion, Shafieloo et al. 10] found that the cosmic expansion acceleration might be slowing 



down, which is different from studies with other SNIa data sets [121 ] . However, once the 
CMB data is included, their result turns out to be consistent with the ACDM model 
very well and the universe is undergoing an accelerating expansion with an increasing 
acceleration. So, there appears some tension between low redshift data (Constitution 
SNIa+BAO) and high redshift (CMB) one. Surprisingly, further analysis using a sub- 
sample (SNLS+ESSENCE+CfA) of the Constitution SNIa reveals that the outcome that 
the cosmic acceleration has been over the peak does not rely on whether the CMB data 
is added, and the tension between SNIa and CMB is reduced si gnifi cantly. Actually, al- 



though previous SNIa data sets, such as Gold06 [13] and Union 



14j . do not support the 



result that the cosmic acceleration is slowing down, the tension between them and CMB 



has already been found [15|, [l6|, and in Ref. |l6|, Nesseris and Perivolaropoulos proposed 



a simple method to find the outliers responsible for it. 

Recently, the largest and latest SNIa sample (Union2) was released by the Supernova 
Cosmology Project (SCP) Collaboration [l^]. It consists of 557 data points. We list the 
subsets in detail in Tab. ()T]). In this Letter, we plan to reexamine the cosmic expanding 
history from the Union2, BAO and CMB by using the popular CPL parametrization. The 
tension between low redshift data and high redshift one is also analyzed in detail. 



II. OBSERVATIONAL DATA 



For SNIa data, we use the latest Union2 compilation released by the Supernova Cos 



mology Project (SCP) Collaboration recently 17J. It consists of 557 data points and is 
the largest published SNIa sample up to now. The statistics of each subset with 3cr outlier 
rejection are detailed in Tab. (JTJ) . We fit the SNIa with cosmological models by minimizing 
the x 2 value of the distance modulus 

557 

X 2 = l^fe) ~ Hobs{Zi)]C^{Zi, Z^lniZj) - Hobs(Zj)}, (2) 

where fj,(z) = 5 log 10 [<iL (z) /Mpc) + 25 is the theoretical value of the distance modulus, 
/i fc s is the corresponding observed one, and C sn (zi, Zj) is the covariance matrix, which was 
detailed in Ref. [l^] and can be found on the web site l . In the present Letter, we will 
use two different covariance matrices, which correspond to the cases with and without 
systematic error, respectively. The luminosity distance cIl(z) is 

For the CPL parametrization, w = wq + wij^, 

E\z) = flo^l + + (1 " Oo™)(l + z)^ 1+w "+^ exp ( - , (4) 



l + z 

where Qom is the present dimensionless density parameter of matter. 



1 http:/ /supernova. lbl.gov/Union/ 
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Since H is a nuisance parameter, we marginalize over it by minimizing the following 
expression 



^ n _ X( > [^.aA'n^^^-lnlO/S]- /In 10 / 2vr \ 

(5) 

to obtain the constraint from SNIa, where oij = /i^ — 25 — 5 log 10 [if ^L(^)]- 

The BAO data considered in our analysis is the distance ratio obtained at z — 0.20 and 
z = 0.35 from the joint analysis of the 2dF Galaxy Redsihft Survey and SDSS data 
which is a relatively model independent quantity and can be expressed as 

D v (z = 0.35) 



D v (z = 0.20) 

with 

2-i 1/3 

D v (z BAO ) 



1.736 ±0.065, (6) 

(7) 



zbao ( f ZBAO 



Performing \ 2 statistics as follows 

2 [D v (z = 0.35)/D v (z = 0.20) - 1.736] 2 

Xbao ~ , (8) 

one can obtain the constraint from BAO. A result from the combination of SNIa and 
BAO is given by calculating XsNia + Xbao- 

Furthermore, in our analysis we add the CMB redshift parameter 
reduce distance at z\ s = 1090 



3, 



which is the 



Zls 

R=V^ -— = 1.71 ±0.019. (9) 



o 



We also apply the x 2 

Y 2 - [it! - L71]2 (10) 
Xcmb ~ Q.019 2 ' 1 1 

to find out the result from CMB and the constraints from SNIa±BAO±CMB are given 
by XsNia ± Xbao + Xcmb- 

III. RESULTS 

We first investigate the constraints on model parameters and then analyze the evo- 
lutionary behavior of the decelerating parameter and Om(z) to probe the properties of 
dark energy and the cosmic expanding history. 



Fig. (OQ) shows the fitting contours of model parameters at the 68.3% and 95% confi- 
dence levels. In the left panel, the systematic error in the SNIa data is ignored, whereas 
in the right panel, it is considered. The dashed, solid and thick solid lines represent the 
results obtained from Union2, Union2+BAO and Union2+BAO+CMB, respectively. The 
point at wo = —1, W\ — denotes the spatially flat ACDM model. We find that, inde- 
pendent of whether the systematic error is taken into account, the outcome from Union2 
is well consistent with that from Union2+BAO, and both Union2 and Union2+BAO ex- 
clude the spatially flat ACDM Universe at 95% confidence level. However, compared to 
the good overlap between regions from Union2 and Union2+BAO, the one obtained from 
Union2+BAO+CMB is relatively isolated and consistent with the ACDM, which means 
that there exists a tension between low redshift data and high redshift one. Obviously, if 
we use the SNIa with the systematic error, this tension is weaker than that from the SNIa 
without. That is, a consideration of systematic errors in the SNIa alleviates this tension 
markedly. 



The evolutionary behaviors of q(z) and Om(z) at the 68.3% confidence level recon- 
structed from Union2+BAO (upper panels) and Union2+BAO+CMB (lower panels) are 
shown in Fig. ([2]). The gray regions and the regions between two long dashed lines repre- 
sent the results without and with the systematic errors in the SNIa, respectively. It is easy 
to see that, for both the SNIa with systematic error and without, there is an apparent 
rise of the values of Om(z) and q(z) in redshifts z < 0.3 for Union2+BAO (upper panels), 
which means that the cosmic acceleration is slowing down. However, this result changes 
dramatically with the addition of CMB in the analysis, as shown in the lower panels of 
Fig. (T5]), which still supports an expansion with an increasing acceleration. These results 



are the same as that derived from the Constitution SN la 



id- 



In order to reduce the tension between low redshift data and high redshift one, Shafieloo 



et al. 



10| use a subsample of Constitution SNIa sample, which is obtained by excluding 



the Gold data, the high z Hubble Space Telescope data and older SNIa data sets in 



-1.5 -1.0 -0.5 0.0 -2.0 -1.5 -1.0 -0.5 0.0 

FIG. 1: The 68.3% and 95% confidence level regions for wq versus w\. In the left panel, the 
system error in the SNIa is ignored, while in the right panel, it is considered. The dashed, 
solid and thick solid lines represent the results obtained from Union2, Union2+BAO and 
Union2+BAO+CMB, respectively. The point at wq = — 1, w\ = represents the spatially 
flat ACDM model. 

Constitution and thus it consists only of SNLS, ESSENCE and CfA. They found that 
the tension is reduced significantly, the outcome does not rely on whether the CMB data 
is added and the cosmic acceleration has been over the peak. Here, we do a similar 



analysis as that in Ref. |l0j by using a subsample of the Union2. This subsample, labeled 
as "Union2S", could be obtained by excluding the Gold data, the high z Hubble Space 
Telescope data and older SNIa data sets in the Union2. It contains 388 data points and 
is given in detail in Tab. §B) (underlined subsets). The fitting contours for wo — w\ and 
reconstructed q(z) and Om(z) are shown in Figs. (l3lkTl) . respectively. From Fig. fl3]), one can 
see that the tension between low redshift data and high redshift one is reduced noticeably, 
and the ACDM is consistent with Union2S with systematic error (Union2S(sys)) and 
Union2S(sys)+BAO+CMB at the 68% confidence level. The left panel of Fig. fl4]) shows 
that, for the case with the systematic error in the SNIa ignored, the evolution of q(z) and 
Om(z) reconstructed using Union2S+BAO is similar to that from Union2S+BAO+CMB. 
Both of them favor that the cosmic acceleration is slowing down. So, the same conclusion 
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FIG. 2: The evolutionary behaviors of q(z) and Om(z) at the 68.3% confidence level. The gray 
regions and the regions between two long dashed lines show the results without and with the 
systematic errors in the SNIa, respectively. The upper and lower panels represent the results 
reconstructed from Union2+BAO and Union2+BAO+CMB, respectively. 



as that from the Constitution SNIa 



lOj is obtained. However, once the systematic error 



in the SNIa is considered, the results from Union2S(sys)+BAO+CMB show that the peak 
of q(z) at z < 0.3 disappears, although Union2S(sys)+BAO still favor a slowing down of 
the present cosmic acceleration. 
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w„ w a 

FIG. 3: The 68.3% and 95% confidence level regions for wq versus w\. A subsample (Union2S) 
of Union2 obtained with the method in [10(] is considered. In the left panel, the system error in 
the SNIa is ignored, while in the right panel, it is considered. The dashed, solid and thick solid 
lines represent the results obtained from Union2S, Union2S+BAO and Union2S+BAO+CMB, 
respectively. The point at wq = — 1, w\ = represents the spatially flat ACDM model. 



Let us now discuss another method in selecting a subsample of the SNIa data, 
which is based upon different considerations. This method was proposed by Nesseris 
and Perivolaropoulos [16j to find the outliers responsible for the tension in the SNIa 
data. In this method, a truncated version of the SNIa can be obtained by calculating 
the relative deviation to the best fit ACDM prediction and adopting a reasonable cut 
\^obs — Hkcdm\/ &obs beyond 1.9 a. Using this method, we find that 39 SNIa points dis- 
tributed in the whole Union2 dataset should be discarded. The names of these 39 SNIa 
are listed in Tab. (JTT1) . Thus, there remain 518 data points and we call them "Union2T". 
The results from Union2T are shown Figs. (0 [6]). We find, from Fig. ([5]), that the ten- 
sion is also reduced significantly, and, the ACDM is consistent with Union2T(sys) and 
Union2T(sys)+BAO+CMB at the 68% confidence level. If the systematic error is ignored, 
the observation favors an expansion with an increasing acceleration at the present once 
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FIG. 4: The evolutionary behaviors of q(z) and Om{z) at the 68.3% confidence level. A sub- 
sample (Union2S) of Union2 obtained with the method in 1(| is considered. The gray regions 
and the regions between two long dashed lines show the results without and with the systematic 
errors in the SNIa, respectively. The upper and lower panels represent the results reconstructed 
from Union2S+BAO and Union2S+BAO+CMB, respectively. 



the CMB is added, although Union2T+BAO still support that the cosmic acceleration 
is slowing down. This can been seen by looking at the grey regions in Fig. ([6]). This 
result is similar with that obtained from Union2, but is different from that from Union2S. 
However, once the systematic errors are considered in the SNIa, both Union2T(sys)+BAO 
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and Union2T(sys)+BAO+CMB favor an expansion with an increasing acceleration at the 
present, which is different from that from Union2 and Union2S. In Tab. (IHIj) . we give the 
X 2 /dof (dof: degree of freedom) value of different datasets, from which, one can see that 
only in the case of Union2T is x 2 /dof significantly improved. That is, according to the 
X 2 /dof criterion, the method proposed in [if]] is preferred. 




-US -1.0 -0.5 0.0 -2.0 -US -1.0 -0.5 0.0 



FIG. 5: The 68.3% and 95% confidence level regions for wo versus w\. A subsample (Union2T) 
of Union2 obtained with the method in [Iff] is considered. In the left panel, the system error in 
the SNIa is ignored, while in the right panel, it is considered. The dashed, solid and thick solid 
lines represent the results obtained from Union2T, Union2T+BAO and Union2T+BAO+CMB, 
respectively. The point at Wq = — 1, w\ = represents the spatially fiat ACDM model. 



IV. CONCLUSION 

In this Letter, we have examined the cosmic expanding history from the latest 558 
Union2 SNIa together with BAO and CMB data. For the SNIa, the data with and 
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FIG. 6: The evolutionary behaviors of q(z) and Om{z) at the 68.3% confidence level. A sub- 
sample (Union2T) of Union2 obtained with the method in [16j is considered. The gray regions 
and the regions between two long dashed lines show the results without and with the systematic 
errors in the SNIa, respectively. The upper and lower panels represent the results reconstructed 
from Union2T+BAO and Union2T+BAO+CMB, respectively. 

without the systematic error are analyzed respectively. The popular CPL parametrization 
is considered. We first find that, independent of whether or not the systematic error 
is considered, there exists a tension between low redshift data (SNIa+BAO) and high 
redshift one (CMB), but for the case with the systematic error considered this tension is 
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weaker than that from the SNIa without. By reconstructing the curves of q(z) and Om(z) 
from Union2+BAO, we obtain that for both the SNIa with and without the systematic 
error the cosmic acceleration has already peaked at redshift z ~ 0.3 and is decreasing. 
However, when the CMB data is added in our analysis, this result changes dramatically 
and the observation favors a cosmic expansion with an increasing acceleration, which 
further confirms the existence of the tension. 

In order to reduce this tension, two different methods given in Refs. Q,Q are con- 
sidered. With the method in [loj], we obtain a subsample of Union2 labeled as Union2S, 
which is given by excluding the Gold data, the high z Hubble Space Telescope data and 
the older SNIa data sets in Union2. Thus 388 data points are left in Union2S. Using 
Union2S, we find that the tension between SNIa+BAO and CMB is reduced markedly. 
For the case without the systematic error both Union2S+BAO and Union2S+BAO+CMB 
favor a decreasing of the cosmic acceleration at z < 0.3. However, once the system- 
atic error is added, Union2S+BAO+CMB support a present increasing cosmic accel- 
eration, although the result from Union2S+BAO is similar with the case without the 
systematic error. According to the method given in If]], we cut 39 data points in 
Union2. Thus, a subsample (Union2T) containing 518 data points is obtained. With 
Union2T the tension is also reduced noticeably. However, when the systematic error is 
ignored, the results from reconstructing q(z) and Om(z) are similar to that given by 
Union2. Union2T+BAO favor that the accelerating expansion of the Universe is slowing 
down, while Union2T+BAO+CMB do not. If the systematic error is considered, both 
Union2T+BAO and Union2T+BAO+CMB support a present increasing cosmic accelera- 
tion. Therefore, when the systematic error in the SNIa is ignored, Union2S and Union2T 
give totally different results once the CMB is added, with one suggesting a slowing-down 
cosmic acceleration, the other just the opposite, although both of them can reduce the 
tension between low redshift data and high redshift one. If the systematic errors in 
the SNIa is considered, we find that the similar results are obtained from Union2S and 
Union2T. Both Union2S+BAO+CMB and Union2T+BAO+CMB support an increasing 
of the present cosmic acceleration. So, in order to have a clear-cut answer on whether the 
cosmic acceleration is slowing down or not, we still need to wait for more consistent data 
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and more reliable methods to analyze them. 

In Tab. flHIl) . the values of x 2 /d°f from different dataset are given. From which one can 
see that this value is significantly improved when Union2T is used. Thus, the XMi n /d°f 
criterion indicates that the method given in [l6] is favored by observations. 

Finally, we must point out that all our results obtained in the present Letter are based 
on the CPL parametrization. If one uses a different parameterization, as in Refs. [lO] , the 
results might change. 
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Set 


Ccwt — 3 


N 


CT sys (68%) 


RMS (68%) 


Hamuy et al. (1996) [19] 


18 


1 5 +0 - 05 

u - lo -0.03 


n i 7 +0.03 
u - 1 ' -0.03 


Krisciunas et al. (2005) [20J 


6 


04+ 13 


11 +0 - 03 

u - li -0.03 


Riess a/. (1999) [21] 


11 


1 5 +0 - 07 
u - lo -0.03 


n i 7 +0.03 
u - 1 '-0.04 


Jha a/. (2006) [22] 


15 


21+ 07 

u ' zi -0.04 


n 99 +0.04 
u - zz -0.04 


Kowalski et al. (2008) [23] 
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07+ - 09 

u - u ' -0.06 


1 5 +0 - 03 
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Hicken et a/. (2009) [7] 
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1 5 +0 - 02 


1Q+ - 01 

u - 1J -o.oi 


Holtzman et al. (2009) [24] 
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in+ - 01 


1 5+ - 01 
u - 10 -o.oi 


Riess et a/. (1998) + HZT [1] 
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31+ 019 


u - oz -0.12 


Perlmutter et a/. (1999) [1] 
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u '^ i -0.09 


64+ - 07 

u - D4 -0.08 


Barris et al. (2004) [25] 
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U - i8 -0.10 


n qo+o.06 

U -3o_0.07 


Amanullah et al. (2008) [26] 
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19+ - 21 
u - ly -0.06 
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u - zi -0.07 


Knop et al. (2003) [27J 
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O 05+ 010 


1 5 +0 - 03 
u - lo -0.02 


Astier et al. (2006) [28] 
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1O+0.03 


0.21+™ 2 


Miknaitis et al. (2007) [29] 
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n 1 q+o- 04 

U ' iy -0.03 


0.29j£g 


Tonry et al. (2003) [30j 
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1 5 +0 - 21 
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Riess et al. (2007) [31] 


31 


1fi+ 006 


u - 4O -0.06 


Amanullah et al. (2010) [17] 


6 


O 00+ 00 
u - uu -0.00 


oo +0 - 00 

u - uu -0.00 


Total 


557 



TABLE I: Statistics of each subset with 3a outlier rejection for Union2 compilation. The 
Union2S consists of the underlined subsets. 
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The subset (39 SNIa) cut in Union2 



1998dx, 1999bm, 2001v, 2002bf, 2002hd, 2002hu, 2002jy, 2003ch, 2003ic, 2006br, 
2006cm, 2006cz, 2007ca, 10106, 200511, 20051p, 2005fp, 2005gs, 2005gr, 2005hv, 
2005ig, 2005iu, 2005jj, 1997k, 2001jm, 1998ba, 03D4au, 04D3cp, 04D3oe, 03D4cx, 
03Dlco, d084, el40, f308, g050, gl20, ml38, 05Str, 2002fx 



TABLE II: The names of SNIa cut in the Union2 with the method in Ref. 



Dataset 


X 2 /dof 


Dataset 


X 2 /dof 


Union2+BAO 


0.962 


Union2(sys)+BAO 


0.938 


Union2+BAO+CMB 


0.964 


Union2(sys)+BAO+CMB 


0.938 


Union2S+BAO 


0.956 


Union2S(sys)+BAO 


0.944 


Union2S+BAO+CMB 


0.955 


Union2S(sys)+BAO+CMB 


0.942 


Union2T+BAO 


0.653 


Union2T(sys)+BAO 


0.631 


Union2T+BAO+CMB 


0.653 


Union2T(sys)+BAO+CMB 


0.630 



TABLE III: Summary of the x 2 /dof from different data sets. 
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